Rapidly detecting disorder in rhythmic biological signals: A spectral entropy measure 

to identify cardiac arrhythmias 
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We consider the use of a running measure of power spectrum disorder to distinguish between the 
normal sinus rhythm of the heart and two forms of cardiac arrhythmia: atrial fibrillation and atrial 
flutter. This spectral entropy measure is motivated by characteristic differences in the power spectra 
of beat timings during the three rhythms. We plot patient data derived from ten-beat windows on 
a "disorder map" and identify rhythm-defining ranges in the level and variance of spectral entropy 
values. Employing the spectral entropy within an automatic arrhythmia detection algorithm enables 
the classification of periods of atrial fibrillation from the time series of patients' beats. When the 
algorithm is set to identify abnormal rhythms within 6 s it agrees with 85.7% of the annotations 
of professional rhythm assessors; for a response time of 30 s this becomes 89.5%, and with 60 s 
it is 90.3%. The algorithm provides a rapid way to detect atrial fibrillation, demonstrating usable 
response times as low as 6 s. Measures of disorder in the frequency domain have practical significance 
in a range of biological signals: the techniques described in this paper have potential application for 
the rapid identification of disorder in other rhythmic signals. 
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I. INTRODUCTION 

Cardiovascular diseases are a group of disorders of the 
heart andblood vessels and are the largest cause of death 
globally 1]. An arrhythmia is a disturbance in the nor- 
mal rhythm of the heart and can be caused by a range of 
cardiovascular diseases. In particular, atrial fibrillation 
is a common arrhythmia affecting 0.4% of the population 
and 5%-10% of those over 60 years old [2] ; it can lead to a 
very high (up to 15-fold) risk of stroke Q. Heart arrhyth- 
mias are thus a clinically significant domain in which to 
apply tools investigating the dynamics of complex bio- 
logical systems Since the pioneering work of Aksel- 
rod et al. on spectral aspects of heart rate variability 
0, such approaches have tended to focus on frequencies 
lower than the breathing rate. By contrast, we develop 
a spectral entropy measure to investigate heart rhythms 
at higher frequencies, similar to the heart rate itself, that 
can be meaningfully applied to short segments of data. 

Conventional physiological measures of disorder, such 
as approximate entropy (ApEn) and sample entropy 
(SampEn) , typically consider long time series as a whole 
and require many data points to give useful results 
0. With current implant technology and the increasing 
availability of portable electrocardiogram (EGG) devices 
, a rapid approach to fibrillation detection is both pos- 
sible and sought after. Though numerous papers propose 
rapid methods for detecting atrial fibrillation using the 
EGG ( 8] and refs therein), less work has been done using 
only the time series of beats or intervals between beats 
[RR intervals). In one study, Tateno and Glass use a sta- 
tistical method comparing standard density histograms 
of Ai?i? intervals [9|. The method requires around 100 
intervals to detect a change in behavior and thus may 
not be a tool suitable for rapid response. 

Measures of disorder in the frequency domain have 



practical significance in a range of biological signals. The 
irregularity of electroencephalography (EEG) measure- 
ments in brain activity, quantified using the entropy of 
the power spectrum, has been suggested to investigate lo- 
calized desynchronization during some mental and motor 
tasks [l^. Thus, the techniques described here have po- 
tential application for the rapid identification of disorder 
in other rhythmic signals. 

In this paper we present a technique for quickly quan- 
tifying disorder in high frequency event series: the spec- 
tral entropy is a measure of disorder applied to the power 
spectrum of periods of time series data. By plotting pa- 
tient data on a disorder map, we observe distinct thresh- 
olds in the level and variance of spectral entropy values 
that distinguish normal sinus rhythm from two arrhyth- 
mias: atrial fibrillation and atrial flutter. We use these 
thresholds in an algorithm designed to automatically de- 
tect the presence of atrial fibrillation in patient data. 
When the algorithm is set to identify abnormal rhythms 
within 6 s it agrees with 85.7% of the annotations of pro- 
fessional rhythm assessors; for a response time of 30 s this 
becomes 89.5%, and with 60 s it is 90.3%. The algorithm 
provides a rapid way to detect fibrillation, demonstrating 
usable response times as low as 6 s and may complement 
other detection techniques. 

The structure of the paper is as follows. Section |ll] in- 
troduces the data analysis and methods employed in the 
arrhythmia detection algorithm, including a description 
of the spectral entropy and disorder map in the context 
of cardiac data. The algorithm itself is presented in Sec. 
HID along with results for a range of detection response 
times. In Sec. IIVI we discuss the results of the algo- 
rithm and sources of error, and compare our method to 
other fibrillation detection techniques. An outline of fur- 
ther work is presented in Sec. |Vl with a summary of our 
conclusions closing the paper in Sec. IVII 



2 



II. DATA ANALYSIS 

After explaining how we symbolize cardiac data in Sec. 
Ill Al the spectral entropy measure is introduced (Sec. 
IIIBI) and appropriate parameters for cardiac data are 
selected (Sec. Ill C[) . We then show how the various 
rhythms of the heart can be identified by their position 
on a disorder map defined by the level and variance of 
spectral entropy values (Sec. IIIPp . 

Data were obtained from the MIT-BIH atrial fibril- 
lation database (afdb), which is part of the physionet 
resource This database contains 299 episodes of 

atrial fibrillation and 13 episodes of atrial flutter across 
25 subjects (henceforth referred to as "patients"), where 
each patient's Holter tape is sampled at 250 Hz for 10 h. 
The onset and end of atrial fibrillation and flutter were 
annotated by trained observers as part of the database. 
The timing of each QRS complex (denoting contraction 
of the ventricles and hence a single, "normal" , beat of the 
heart) had previously been determined by an automatic 
detector [12| . 



A. Symbolizing cardiac data 

We convert event data into a binary string, a form 
appropriate for use in the spectral entropy measure. 
The beat data is an event series: a sequence of pairs 
denoting the time of a beat event and its type. We 
categorize normal beats as N and discretize time into 
short intervals of length r (for future reference, sym- 
bols are collected with summarizing descriptions in Table 
m. Each interval is categorized as or depending on 
whether it contains no recorded event or a normal beat, 
respectively. This yields a symbolic string of the form 
...000N00N0N000N.... This symbolic string can 
be mapped to a binary sequence (A^ ^ 1, ^ 0). This 
procedure is shown schematically in Fig. [T] Naturally, 
this categorization can be extended to more than two 
states and applied to other systems. For example, ec- 
topic beats (premature ventricular contractions) could be 
represented by V to yield a symbolic string drawn from 
the set {0,N, V}. An additional map could then be used 
to extract a binary string representing the dynamics of 
ectopic beats. 



B. Spectral entropy 

We now present a physiological motivation for using a 
measure of disorder in the context of cardiac dynamics, 
followed by a description of the spectral entropy mea- 
sure. Following Ref. 0], atrial fibrillation is character- 
ized by the physiological process of concealed conduction 
in which the initial regular electrical impulses from the 
atria (upper chamber of the heart) are conducted inter- 
mittently by the atrioventricular node to the ventricles 
(lower chamber of the heart). This process is responsi- 



ble for the irregular ventricular rhythm that is observed. 
Atrial flutter has similar causes to atrial fibrillation but 
is less common; incidences of flutter can degenerate into 
periods of fibrillation. Commonly, alternate electrical 
waves are conducted to the ventricles, maintaining the 
initial regular impulses originating from the atria. This 
results in a rhythm with pronounced regularity. Nor- 
mal sinus rhythm can be characterized by a slightly less 
regular beating pattern occurring at a slower rate com- 
pared to atrial flutter. Example electrocardiograms for 
the three rhythms are shown in the boxed-out areas of 
Fig. [3l below. 

Given these physiological phenomena, the spectral en- 
tropy can be used as a natural measure of disorder, en- 
abling one to distinguish between these three rhythms of 
the heart. Presented with a possibly very short period 
of heart activity one can create a length-L, duration-Lr, 
binary string. One then obtains the corresponding power 
spectrum of this period of heart activity using the dis- 
crete Fourier transform |13i] . Given a (discrete) power 
spectrum with the ith frequency having power Q, one 
can define the "probability" of having power at this fre- 
quency as 



When employing the discrete Fourier transform, the sum- 
mation runs from i = 1 to i = ^. One can then find the 
entropy of this probability distribution [with i having the 
same summation limits as in Eq. ([1])] : 

H = ^-p^\og2P^. (2) 

i 

Breaking the time series into many such blocks of dura- 
tion Lr, each with its own spectral entropy, thus returns 
a time series of spectral entropies. Note that this measure 
is not cardiac specific and can be applied to any event 
series. For example, a sine wave having period an integer 
fraction of Lt will be represented in Fourier space by a 
delta function (for Lt — > oo) centered at its fundamental 
frequency; this gives the minimal value for the spectral 
entropy of zero. Other similar frequency profiles, with 
power located at very specific frequencies, will lead to 
correspondingly low values for the spectral entropy. By 
contrast, a true white noise signal will have power at all 
frequencies, leading to a flat power spectrum. This case 
results in the maximum value for the spectral entropy: 




As will be discussed in the following section, H can be 
normalized by Hmax to give spectral entropy values in 
the range [0,1]. 

Note that analytical tools relying on various interbeat 
intervals have been devised in the past (e.g. d, [13, Elj)- 
Here, we demonstrate how our measure relates to those 
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studies. Any series of events can be represented by 

/(i) = E<5(t-t,), (4) 

k 

where tk is the time when an event (beat) occurs. The 
corresponding power spectrum is, then, 

P{Lj)^Y.cos{uj\tk-tk'\). (5) 

k,k' 

The spectral entropy is, by definition, 

Hcont. ^ J dujp{uj) log p{lj), (6) 

where p{uj) = Pii^)/ J du!'P{LLj'). We therefore see that 
Eq. ([6]) depends on all of the intervals between any two 
events (c.f. Eq. ([5])). This is in contrast to studies 
on the distribution of beat-ncxt-beat intervals in p^ . 
We believe that this generalization enriches the struc- 
ture captured in the short-time segments and thus al- 
lows for the shortening of the detection response time in 
our arrhythmia detection algorithm. We finally note that 
since the spectral entropy depends only on the shape of 
the power spectrum, it is relatively insensitive to small, 
global, shifts in the spectrum of the signal. 

C. Parameter selection 

We now introduce parameters for the spectral entropy 
measure and select values appropriate for cardiac data. 
The sampling interval acts like a low pass-filter on the 
data since all details at frequencies above 1/(2t) Hz, the 
upper frequency limit, are discarded [l^. The sampling 
interval must be sufficiently small such that no useful 
high-frequency components are lost. We choose a sam- 
pling interval r=30 ms, since processes like the heart beat 
interval, breathing and blood pressure fluctuations occur 
at much lower frequencies. The upper frequency limit in 
the power spectrum is consistent with the inclusion of all 
dominant and subsidiary frequency peaks present during 
atrial fibrillation [13]. 

We call the duration over which the power spectrum 
is found, and hence a single spectral entropy value is ob- 
tained, the spectral entropy window: a — Lt {L is the 
number of sampling intervals required). With our value 
for T, the shortest spectral entropy window giving suffi- 
cient resolution in the frequency domain for cardiac data 
is found for L around 200, a ~ 6 s. This value for a is 
equivalent to approximately ten beats on average over the 
entire afdb. It is consistent with previous work on animal 
hearts looking at the minimum window length required 
to determine values for the dominant frequencies present 
during atrial fibrillation To take into account the 

heterogeneity of patients' resting heart rates (HRs), we 
fix T and use an L value for each patient such that there 
are on average 10 beats in each spectral entropy window. 



Thus, a = L{HR)t — a{HR). All subsequent parame- 
ters that are determined by L will similarly be a function 
of the average heart rate; we will henceforth drop the 
HR notation for clarity, with the dependence on average 
heart rate understood implicitly. Patients with higher 
average heart rate require smaller L, and therefore have 
a shorter spectral entropy window. By invoking individ- 
ual values for L, the maximum spectral entropy for each 
patient is constrained to a particular value: Hmax [c.f. 
Eq. ([3])]. To make spectral entropy values comparable, 
we normalize the basic spectral entropy values for each 
patient [Eq. ^] by their theoretically maximal spectral 
entropy value. The spectral entropy can thus take values 
in the range [0,1]. In choosing L near its minimally us- 
able value, we necessarily have a small number of beats 
compared to the window length a. In such cases, a win- 
dow shape having a low value for the equivalent noise 
bandwidth (ENBW) is preferable [H]. The ENBW is a 
measure of the noise associated with a particular window 
shape: it is defined as the width of a fictitious rectan- 
gular filter such that power in that rectangular band is 
equal to the actual power of the signal. The condition 
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FIG. 1: Schematic of cardiac data analysis and the auto- 
matic arrhythmia detection algorithm. A full description of 
the data analysis (stages 1-3) is given in the Data Analysis 
section, Sec. [Ill of the text; the remaining steps (stages 4-5) 
are described in the Algorithm section, Sec. IIIII MIT-BIH 
atrial fibrillation database event data (stage 1) are discretized 
at sampling interval r, then mapped to give a binary series 
representing the dynamics of regular beats A'' (stage 2). A 
sequence of spectral entropy windows, of length a, is applied 
with overlap parameter a to obtain a series of spectral entropy 
values (stage 3). Variance windows, of length /3 with overlap 
parameter b, are applied to obtain a series of variance values. 
Threshold conditions in the level and variance of spectral en- 
tropy values allows for the classification of periods of atrial 
fibrillation (AF) and other rhythms (N), typically normal si- 
nus rhythm (stage 4). Finally, the most frequent prediction 
in each modal smoothing window, of length 7 with overlap 
parameter c, is identified {AF', N'} to obtain the final algo- 
rithm output (stage 5). For definitions and typical values for 
algorithm parameters, see Table IH 
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Window 


Symbol 


Definition 


Typical value 


Overlap 


Typical value 


Spectral entropy 


a 


Lt 


6 s 


a 


1.5 s 


Variance 


P 


Ma = MLt/A 


6 s, 30 s, 60 s 


b 


1.5 s 


Modal smoothing 


7 


2P + h= (2M + l)Lr/4 


12 s, 60 s, 120 s 


c 


1.5 s 



TABLE I: Summary of arrhythmia detection algorithm window and overlap symbols. A full description of the spectral entropy 
and variance windows is given in the Data Analysis section, Sec. UT] of the text; the modal smoothing window is described in 
the Algorithm section, Sec. IIIII Cardiac data in the MIT-BIH atrial fibrillation database is sampled at intervals of t—30 ms. 
The number of intervals contained in the spectral entropy window, L, is chosen for each patient such that the spectral entropy 
window is expected to contain ten beats. In the variance window, M represents the number of spectral entropy values used in 
finding the variance; for response times 6 s, 30 s, 60 s, we consider M equal to 4, 20, 40, respectively. Specifying r, L and M 
fixes the remaining parameters. We define overlap parameter a=Q/4. For simplicity, we set c=b=a. 



for low ENBW is satisfied by the rectangular window. To 
maximize the available data, a sequence of overlapping 
rectangular windows separated by a time a is used. This 
results in a series of spectral entropy values also sepa- 
rated by a. We follow the convention of using adjacent 
window overlap of 75% [l9| , leading to a window separa- 
tion time: a = Lt/A. This gives a typical value for a of 
1.5 s. A summary of window and overlap parameters is 
presented in Table [H 

Figure [5] illustrates the spectral entropy measure ap- 
plied to patient 08378 from the afdb. We identify three 
distinct levels in the spectral entropy value correspond- 
ing to the three rhythms of the heart assessed in the 
annotations. Beating with a relatively regular pattern, 
which can be associated with normal sinus rhythm, sets 
a baseline for the spectral entropy. The irregularity as- 



sociated with fibrillation causes an increase in the value, 
with the pronounced regularity of flutter identifiable as 
a decrease in the spectral entropy. We note that power 
spectrum profiles in frequency space should remain qual- 
itatively similar for a given rhythm type, regardless of 
the underlying heart rate. For example, periodic signals 
can be characterized by peaks at constituent frequen- 
cies, independent of the beat production rate; similarly, 
highly disordered signals can be consistently identifiable 
by their flat power spectra. This confers a significant ad- 
vantage over methods relying solely on the heart rate. We 
find considering only the instantaneous heart rate and its 
derivatives to be insufficient in consistently distinguish- 
ing between sinus rhythm, fibrillation and flutter; this 
point is addressed further in the Discussion section (Sec. 
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FIG. 2: (Color online) Spectral entropy time series (top), 
professional rhythm annotation (middle) , and arrhythmia de- 
tection algorithm prediction (bottom) for patient 08378 from 
the MIT-BIH atrial fibrillation database. Event data is sam- 
pled at 30-ms intervals approximately 200 times such that 
there are on average ten beats per spectral entropy window. 
Each window, of length 6 s for a typical patient, contributes 
one value of the spectral entropy; windows have a typical 
overlap of 1.5 s. For the rhythm annotation and algorithm 
prediction: AF denotes atrial fibrillation, AFL denotes atrial 
flutter, and N represents all other rhythms. The algorithm 
prediction (primed symbols omitted for clarity) demonstrates 
good agreement with professional annotations; shown for a 
response time of 30 s, thresholds: r/i6=0.84, r/i=0.70 and 
-l'/,6=0.018. 



D. Cardiac disorder map 

Having identified differences in the level of the spec- 
tral entropy measure corresponding to different rhythms 
of the heart, we suggest that there should be a similar 
distinction in the variance of a series of spectral entropy 
values. We propose that the fibrillating state may repre- 
sent an upper limit to the spectral entropy measure; once 
this state is reached, variations in the measure's value are 
unlikely until a new rhythm is established. By contrast, 
the beating pattern of normal sinus rhythm is not as dis- 
ordered as possible and can therefore show variation in 
the spectral entropy values taken. Inspecting the data, 
one frequently observes transitions between periods of 
very regular and more irregular (though still clearly si- 
nus) beating. Thus, normal sinus rhythm will naturally 
explore more of the spectral entropy value range than 
atrial fibrillation, which is consistently irregular in char- 
acter (including some dominant frequencies [iTj). Fur- 
thermore, in defining the spectral entropy window to be 
constant for a given patient, some dependence on the 
heart rate is retained, despite accounting for each pa- 
tient's average heart rate. This dependence can cause ad- 
ditional harmonics to appear in the power spectrum, in- 
creasing the variation of spectral entropy values explored 
during normal sinus rhythm. Last, windows straddling 
transitional periods of the heart rate will also demon- 
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FIG. 3: (Color online) Cardiac disorder map for all 25 patients in the MIT-BIH atrial fibrillation database (afdb). Boxed- 
out area: example electrocardiograms for normal sinus rhythm, atrial fibrillation and atrial flutter, taken from patient 04936. 
Spectral entropy values are obtained from windows of event data expected to contain ten beats; data is sampled at 30-ms 
intervals. For a typical patient, each spectral entropy window is around 6 s in length and has an overlap with adjacent windows 
of 1.5 s. For each point on the disorder map, the standard deviation and average spectral entropy level is calculated from M 
adjacent spectral entropy values: we call this the variance window, /3. Here, we have M equal to 20 and so /3 has a typical length 
of 30 s; /? represents the response time of the algorithm. Normalized frequency histograms are disorder map projections onto 
the relevant axes. Rhythm assessments, {N, AF, AFL}, are provided in the afdb. Atrial fibrillation is situated in the upper left 
of the disorder map, consistent with having a high value for the spectral entropy and a low value for the variance. Atrial flutter 
has a lower average value for the spectral entropy, as expected. Fibrillation thresholds for the arrhythmia detection algorithm 
are set at r/ii,=0.84 for the spectral entropy level and <l>/ib=0.018 for the standard deviation, as indicated on the disorder map. 



strate atypical power spectra, further compounding the 
increase in the variance when comparing normal sinus 
rhythm to atrial fibrillation. We do not conjecture on 
(and do not observe) a characteristic difference in the 
variance of spectral entropy values for atrial flutter, re- 
lying on the spectral entropy level to distinguish the ar- 
rhythmia from fibrillation and normal sinus rhythm. 

In theory, the spectral entropy can take values in the 
range [0,1]. Possible variances in sequences of spectral 
entropy values then lie in the range [0,i]. These two 
ranges determine the two-dimensional cardiac disorder 
map. In practice, we plot the standard deviation rather 
the variance for clarity, and so rhythm thresholds are 
given in terms of the standard deviation. Due to finite 
time and windowing considerations, the spectral entropy 
is restricted to a subset of values within its possible range. 
We attempt to find limits in the values that the spectral 
entropy can take by applying the measure to synthetic 



event series: a periodic series with constant interbeat in- 
terval, and a random series drawn from a Poisson prob- 
ability distribution with a mean firing rate. For a heart 
rate range of 50 beats per minute (bpm) to 200 bpm in 
1-bpm increments we obtain 150 synthetic time series for 
the periodic and Poisson cases, respectively. The average 
spectral entropy value over the 150 time series in the pe- 
riodic case is 0.67±0.04; the average value in the Poisson 
case is 0.90±0.01. By assuming the maximum variance 
to occur in a rhythm that randomly changes between the 
periodic and Poisson cases with equal probability, an ap- 
proximate upper bound for the standard deviation can be 
calculated: using the two average spectral entropy values 
in the definition of the standard deviation, we find the 
upper bound to be approximately 0.115. 

Figure [3] illustrates the cardiac disorder map for all 25 
patients comprising the afdb. The standard deviation 
is calculated from M adjacent spectral entropy values 
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(separated by a), corresponding to a duration of /3 = 
Ma — MLt/A; we call /3 the variance window. In this 
case, we have M equal to 20 and so (3 has a length of 
30 s for a typical patient. We will see in the following 
section that (3 sets the response time of the arrhythmia 
detection algorithm. The smallest useable number for 
M is 4, corresponding to the rapid response case where 
f3 is typically 6 s. We have M equal to 40 for the case 
where /3 is typically 60 s. In Fig. [31 each value of the 
standard deviation is plotted against the average value 
of the spectral entropy within the variance window, and 
is colored according to the rhythm assessment provided 
in the annotations. As with spectral entropy windows, 
variance windows have an overlap, b. For simplicity, we 
set b=a, giving a typical value of 1.5 s. Note that b can 
take any integer multiple of a, though doing so does not 
alter the results substantially. 

One observes atrial fibrillation to be situated in the 
upper left of the disorder map, consistent with having a 
high value for the spectral entropy and a low value for the 
variance. Atrial flutter has a lower average value for the 
spectral entropy, as expected. For the given case with (3 
typically 30 s, we determine fibrillation to exhibit spec- 
tral entropy levels above r/ib=0.84, with flutter present 
below r/;=0.70. A standard deviation threshold can be 
inferred at around ^ fib— 0-018, with the majority of flb- 
rillating points falling below that value. Although be- 
yond the expository purpose of this paper, we note that 
these approximate thresholds can be further optimized 
using, for example, discriminant analysis (20| . Disorder 
maps for the three detection response times (6 s, 30 s, 
60 s) are qualitatively similar; increasing the length of 
the variance window improves the separation of rhythms 
in the disorder map at a cost of requiring more data per 
point. From these observations, we hypothesize threshold 
values in the spectral entropy level and variance that dis- 
tinguish the two arrhythmias from normal sinus rhythm. 
In the following section, thresholds drawn from the disor- 
der map are used in an arrhythmia detection algorithm. 



III. ALGORITHM 

In this section, we present a description of the auto- 
matic arrhythmia detection algorithm (Sec. IIII Ap . fol- 
lowed by results for a range of detection response times 
(Sec. MM- 



A. Arrhythmia detection algorithm 

The arrhythmia detection algorithm uses thresholds in 
the level and variance of spectral entropy values observed 
in the cardiac disorder map to automatically detect and 
label rhythms in patient event series data. The afdb con- 
tains significantly fewer periods of atrial fiutter compared 
to atrial fibrillation and normal sinus rhythm (periods of 
fiutter total 1.27 h, whereas periods of fibrillation total 



91.59 h), the typical length of periods of flutter is of the 
order tens of seconds. Of the eight patients annotated as 
having flutter, only patients 04936 and 08378 have pe- 
riods of flutter long enough (i.e., > (3) for analysis by 
the algorithm. For this reason we do not include here 
the fiutter prediction method of the algorithm, although 
extensions including fiutter follow a similar principle and 
are simple in practice to implement. Other studies using 
the afdb (e.g., |9j) restrict themselves to methods dif- 
ferentiating only between fibrillation and normal sinus 
rhythm. Additional comments on the practicality of de- 
tecting atrial fiutter and selected results for fiutter will 
be given in the Discussion section (Sec. IIV A|l . 

The five stages of the algorithm are shown in Fig. [TJ 
The first three stages have been covered in depth as part 
of the Data Analysis section, but we include a brief sum- 
mary here for completeness. We first obtain a binary 
string representing the dynamics of the heart for a given 
patient by discretizing the physionet data every t=30 ms 
(stage 1 to stage 2). In stage 3, the spectral entropy mea- 
sure is applied for windows of duration a = Lt, with L 
chosen for each patient such that there are on average 
ten beats within the spectral entropy window, giving a 
as 6 s for a typical patient. Using an overlap parameter 
a (typically 1.5 s), leads to a series of spectral entropy 
values separated in time by this amount. Given no prior 
knowledge of the provided rhythm assessments, we cal- 
culate the standard deviation and average magnitude of 
M spectral entropy values in variance windows of length 
(3 = Ma preceding a given time point. We use the exam- 
ple case of M equal to 20 (giving /3 as 30 s for a typical 
patient). The level and standard deviation thresholds 
for atrial fibrillation are set consistent with values ob- 
tained from the cardiac disorder map, for this case we 
determine r/ib=0.84 and <I>/,;f,=:0.018. Stage 4 generates 
preliminary predictions for the rhythm state of the heart: 
we denote as fibrillating (AF) instances where the spec- 
tral entropy level is greater than Tfih and the standard 
deviation is less than ^/ib, with all other combinations 
considered to be normal sinus rhythm (N) [2l[. Setting 
the overlap of variance windows such that & = a, we ob- 
tain a string of rhythm predictions drawn from the set 
{AF, N} and separated in time by b. 

Finally, in stage 5 we apply a rudimentary smoothing 
procedure to the initial string of rhythm predictions. For 
a particular prediction, we consider a preceding period 
7 = 2/3-1-6= (2M + 1)Lt/A, leading in this example to a 
typical length for 7 of 61.5 s. We find the modal predic- 
tion: the prediction {AF, N} occurring most frequently 
in 7, labeling the modal prediction {AF', N'}. We call 
7 the modal smoothing window. In this form, we under- 
stand the windows (3 and 7 as setting the response time 
of the algorithm: (3 is defined in terms of the number 
of preceding spectral entropy values required for a given 
prediction; for 7 to register a change in rhythm, over half 
of the predictions must suggest the new rhythm. The re- 
sponse time is then ^, which is approximately equal to 
(3. We have the modal smoothing windows overlapping 
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with parameter c=b=a. This results in a final time series 
of predictions and constitutes the output of the arrhyth- 
mia detection algorithm for a given patient. An example 
of the algorithm output for patient 08378 (including a 
threshold for atrial flutter) is shown in Fig. [2l 

We apply the above steps, comprising the three data 
windows (a, /?, 7), to each patient in the afdb. Specifying 
r, L and M fixes the remaining parameters, their exact 
magnitude determined by L. A summary of windowing 
symbols can be found in Table ID Values for the atrial 
fibrillation threshold parameters (T fib and ^fib) are kept 
the same for each patient for a given response time. The 
results obtained from the algorithm are described in the 
following section. 



B. Algorithm results 



We now present the results of the cardiac arrhythmia 
detection algorithm for atrial fibrillation. The follow- 
ing window parameters were used: r is set to 30 ms, L 
is chosen such that a is expected to contain 10 beats, 
and M is set to 20, windows have overlap parameters 
c=b=a=a/A (for typical patients in the afdb, a ~ 6s, 
P ~ 30s, 7 ~ 61.5s, and a ^ 1.5s). Threshold values 
for fibrillation are set at T fib— 0.8A for the spectral en- 
tropy level and ^ fn,— 0.018 for the standard deviation. 
Each prediction produced by the algorithm (denoted by 
a primed symbol) is compared with the rhythm assess- 
ment documented in the database and can be classified 
into one of four categories H^]: true positive (TP), AF 
is classified as AF'; true negative (TN), non-AF is clas- 
sified as non-AF'; false negative (FN), AF is classified as 
non-AF'; false positive (FP), non-AF is classified as AF'. 
Percentages of these quantities for each patient and for 
the entire afdb are given in Table [TTl Overall, we obtain 
a predictive capability (assessed using the percentage of 
predictions agreeing with the provided annotations) of 
89.5%. The sensitivity and specificity metrics are de- 
fined by TP/(TP+FN) and TN/(TN-hFP), respectively. 
The predictive value of a positive test (PV+) and the 
predictive value of a negative test (PV_) are defined by 
TP/(TP-hFP) and TN/(TN-FFN), respectively. These, 
and results for other values of /3 are given in Table [TTTl 

In repeating the algorithm with different values for 
the variance window, shorter (3 represents a quicker re- 
sponse time. We obtain for each /3 a new disorder map 
to determine the relevant threshold values. For the rapid 
response case, (3 typically 6 s, we alter the fibrillating 
thresholds in the arrhythmia detection algorithm to be 
r/if,=0.855 and <I>jib=0.016; we find a predictive capa- 
bihty of 85.7%. With (3 typically 60 s, the fibrillating 
thresholds become Tfib=0.8A and <i>/ib=0.019; the pre- 



Patient 


TP (%) TN (%) 


FN (%) FP (%) 


00735 


0.8 


85.0 


0.0 


14.2 


03665 


29.8 


30.4 


37.8 


2.0 


04015 


0.5 


92.4 


0.2 


6.9 


04043 


8.9 


76.5 


13.1 


1.5 


04048 


0.4 


98.8 


0.7 


0.1 


04126 


3.3 


78.3 


0.6 


17.8 


04746 


53.6 


43.8 


0.8 


1.8 


04908 


7.0 


88.2 


1.6 


3.2 


04936 


43.1 


19.0 


36.3 


1.6 


05091 


0.0 


85.6 


0.2 


14.2 


05121 


56.9 


30.5 


8.4 


4.2 


05621 


0.9 


94.9 


0.4 


3.8 


06426 


92.7 


1.9 


3.1 


2.3 


06453 


0.4 


97.7 


0.7 


1.2 


06995 


42.8 


47.1 


3.0 


7.1 


07162 


100.0 


0.0 


0.0 


0.0 


07859 


83.1 


0.0 


16.9 


0.0 


07879 


53.3 


38.1 


8.5 


0.1 


07910 


13.5 


85.7 


0.5 


0.3 


08215 


80.0 


19.7 


0.3 


0.0 


08219 


18.3 


59.8 


3.8 


18.1 


08378 


20.0 


77.3 


2.4 


0.3 


08405 


68.9 


28.4 


2.7 


0.0 


08434 


3.8 


91.6 


0.2 


4.4 


08455 


65.6 


31.5 


2.9 


0.0 


Total 


36.1 


53.4 


6.5 


4.0 




True: 


89.5% 


False: 


10.5% 



TABLE II: Results of the arrhythmia detection algorithm us- 
ing data in the MIT-BIH atrial fibrillation database. For the 
parameters used, see Algorithm results section (Sec. IIII B[l . 
Algorithm predictions (primed symbols) are compared to an- 
notated rhythm assessments. TP, AF is classified as AF'; 
TN, non-AF is classified as non-AF'; FN, AF is classified as 
non-AF'; FP, non-AF is classified as AF'. 



M P 



4 6s 
20 30s 
40 60s 



True (%) 



85.7 
89.5 
90.3 



Sens. (%) Spec. (%) PV+ (%) PV- (%) 



82.1 


88.4 


83.9 


87.0 


84.8 


92.9 


89.8 


89.2 


83.6 


95.2 


92.8 


88.7 



fib 

dictive capability is 90.3%. 



TABLE III: Summary of results for variance windows of dif- 
ferent lengths. Length is set by parameter M=4, 20, 40, 
giving durations for typical patients: (3 ~ 6s, 30s, 60s, 
respectively. Shorter (3 indicates a quicker response time. 
Metrics defined as, sensitivity, TP/(TP-|-FN); specificity, 
TN/(TN-HFP); PV+, TP/(TP+FP); PV_, TN/(TN-|-FN). 



IV. DISCUSSION 

We begin with an exposition of the results presented 
in the previous section and the effects of different param- 
eter values on the output of the arrhythmia detection al- 
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gorithm. This is followed by a discussion, with reference 
to the electrocardiograms provided as part of the afdb, 
of disagreements between the provided rhythm annota- 
tions, measures relying solely on the heart rate, and the 
predictions of our algorithm (Sec. IIV Ap . Having shown 
that some of the annotations may be unreliable, we com- 
ment on situations where the algorithm may still present 
incorrect predictions (Sec. IIV Bp . The benefits of the 
spectral entropy measure compared to other fibrillation 
detection methods is then given (Sec. IIVC|) . We close 
the section with a discussion of the systematic windowing 
errors present in our procedure (Sec. lIVPp . 

Instances of atrial fibrillation constitute approximately 
40% of the afdb. If we consider a null-model where we 
constantly predict normal sinus rhythm, we would expect 
a predictive capability of around 60%. In Table [1111 we 
observe an improvement in the predictive capability of 
the detection algorithm when the length of the variance 
window, (3, is increased from 6 s (85.7%) to 60 s (90.3%) 
for a typical patient. The choice of shorter /3 improves 
the response time of the algorithm by requiring less data 
per prediction; values for (3 less than 6 s do not incorpo- 
rate enough data to give meaningful results. Increasing f3 
beyond 30 s improves the predictive capability very little. 
This suggests that additional factors, independent of the 
specific parameters chosen here, need to be considered. 
Results in Table HIl for the case f3 typically 30 s indicates 
an overall predictive capability of 89.5%. For individ- 
ual patients, the predictive capability ranges from 60.2% 
(patient 03665) to 100% (patient 07162). To explain this 
variation, we investigate the form of patient ECGs during 
periods of disagreement between annotation and predic- 
tion. Examples of the ECGs referred to in Sees. IIV Al 
and IIVBI are included in the supplementary information 
that accompanies this paper [23|. 



A. Disagreements with annotations 

Rhythm assessments have been questioned before @; 
here, we give explicit examples where we believe the 
ECGs to suggest a rhythm different from that given by 
the annotation. We observe in the ECGs of patients 
08219 and 08434 periods of atrial fibrillation that we be- 
lieve to have been missed in the annotations but are cor- 
rectly identified by our detection algorithm [2I]. Cases 
such as these serve to negatively impact the results of 
the algorithm unfairly; however, we note that such in- 
stances comprise a small proportion of the afdb. Atrial 
flutter may have been misannotated in patients 04936 
and 08219 [1^; in particular, two considerable periods of 
flutter may have been annotated incorrectly in patient 
04936. This unreliability of rhythm assessment, com- 
pounded with the limited number of periods of atrial 
flutter in the database, prevents us from drawing mean- 
ingful quantitative conclusions regarding the success of 
the detection algorithm in identifying flutter. Despite 
this, we believe that the spectral entropy is in principle 



still capable of identifying flutter (see Fig. [2|) . Returning 
to the two patients with significant periods of flutter, we 
run the algorithm with the inclusion of a threshold for 
atrial flutter motivated by each patient's individual dis- 
order map: Tfi (other parameters as per the Algorithm 
results section with M — 20). For patient 08378 with 
T fi — 0.70, we find 86.3% agreement with the annota- 
tions for flutter; for patient 04936 with T fi = 0.81, we 
find 66.9% agreement, bearing in mind the points raised 
above. 

Consideration of ECGs demonstrates the inability of 
measures relying solely on the heart rate and its deriva- 
tives to consistently distinguish between fibrillation, flut- 
ter and other rhythms. Atrial fibrillation is characteris- 
tically associated with an elevated heart rate (100-200 
bpm) ^] ; atrial flutter exhibits an even higher heart rate 
(>150 bpm) with a sharp transition from normal sinus 
rhythm. This expected behavior, whilst found to hold 
qualitatively for the majority of patients, fails during 
large periods for patient 06453 and is completely reversed 
for patient 08215 f2^. The resting heart rate is also found 
to differ dramatically between patients in the afdb. The 
spectral entropy, being less susceptible to variations in 
the heart rate [24I , is better suited to form the basis of a 
detection algorithm compared to a measure relying solely 
on heart rate. 



B. Other rhythms 

The unreliability of parts of the annotations still does 
not account for all false predictions produced by the 
detection algorithm. We suggest the presence of other 
rhythms within the afdb to be an additional factor that 
needs to be considered. Table IIIII shows the sensitivity 
metric to be consistently lower for all values of f3, suggest- 
ing a bias towards false negatives (FNs occur when AF is 
classified as non-AF'). FNs total 6.5% for /3 typically 30 s 
in Table HH and comprise 36.3% of predictions for patient 
04936. Given our requirement in the detection algorithm 
for periods that are classed as AF to satisfy both a spec- 
tral entropy level and variance condition, FNs are most 
likely to arise when one threshold condition fails to be 
met. Cases where the variance threshold is not satisfied 
may be associated with the physiological phenomena of 
fib- flutter and paroxysmal atrial fibrillation, and would be 
located right of the standard deviation threshold on the 
disorder map (Fig. [3]) . Fib-flutter corresponds to periods 
where the rhythm transitions in quick succession between 
atrial fibrillation and flutter ^28] , with paroxysmal fibril- 
lation describing periods where atrial fibrillation stops 
and starts with high frequency. Such behavior naturally 
causes the variance to increase and one might question 
whether it is still appropriate to classify those periods as 
standard atrial fibrillation. We identify in the EGG of 
patient 04936 periods of fib-flutter which likely accounts 
for the high proportion of FN results [1^; by inspect- 
ing the patient's disorder map, we indeed observe points 
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annotated as atrial fibrillation with uncharacteristically 
high standard deviation, signifying that fib-flutter would 
be a more accurate rhythm classification. Cases where 
the spectral entropy level threshold is not met can occur 
when QRS complexes indicative of atrial fibrillation ap- 
pear with unusually regular rhythm; such behavior would 
lie below the level threshold on the disorder map. Ow- 
ing to the small number of beats contained within each 
window, such occurrences inevitably arise; the process of 
modal smoothing lessens the impact of this phenomenon 
in the arrhythmia detection algorithm. 

False positives (FPs occur when non-AF is classified as 
AF'), which comprise 4.0% of the afdb for f3 typically 30 
s, may also have a physiological explanation. During si- 
nus arrhythmia, there are alternating periods of slowing 
and increasing node firing rate, while still retaining QRS 
complexes indicative of normal sinus rhythm. These al- 
ternating periods increase the irregularity of beats within 
the spectral entropy window. If the variance threshold is 
also satisfied, sinus arrhythmia may be incorrectly classi- 
fied as AF' by the arrhythmia detection algorithm. Sinus 
arrest occurs when the sinoatrial node fails to fire and 
results in behavior that is similar in principle to sinus ar- 
rhythmia; these two conditions are likely responsible for 
the high proportion of FPs (14.2%) that are observed in 
patient 05091 [s^. 



C. Comparison to other methods 

Vikman et al. showed that decreased ApEn values of 
heart beat fluctuations have been found to precede (at 
timescales of the order an hour) spontaneous episodes of 
atrial fibrillation in patients without structural heart dis- 
ease ^l| . We stress that the algorithm presented here is 
not intended to predict in advance occurrences of fibril- 
lation; rather, it is designed to detect the onset of fib- 
rillation as quickly as possible using only interbeat in- 
tervals. Tateno and Glass ^] present an atrial fibrilla- 
tion detection method that is statistical in principle and 
based upon an observed difference in the standard den- 
sity histograms of /S.RR intervals (the difference in suc- 
cessive interbeat intervals) . A series of reference standard 
density histograms characteristic of atrial fibrillation (as 
assessed in the annotations) are first obtained from the 
afdb. Their detection algorithm is re-run on the afdb by 
taking 100 interbeat intervals and comparing them to the 
reference histograms, where appropriate predictions can 
then be made. The reference histograms rely on the cor- 
rectness of the annotations in order to determine fibril- 
lation, whereas the thresholds in our algorithm are only 
weakly dependent on the data set under consideration. 
Figure [3] is an empirical observation, in future analyses 
we would like to use fibrillation thresholds derived from 
a data set separate from the one under consideration. 

Sarkar et al. have developed a detector of atrial fibril- 
lation and tachycardia that uses a Lorentz plot of Ai?i? 
intervals to differentiate between rhythms [32l |. The de- 



tector is shown to perform better for episodes of fibril- 
lation greater than 3 min and has a minimum response 
time of 2 min. By contrast, our method is applicable to 
short sections of data, enabling quicker response times 
to be used. We see our algorithm complementing other 
detection techniques, with the potential for an implemen- 
tation that combines more than one method. Combin- 
ing methods becomes increasingly relevant when running 
algorithms on data sets containing a variety of arrhyth- 
mias. As noted in Q, other arrhythmias often show ir- 
regular RR intervals, and previous studies have found 
difficulty in detecting atrial fibrillation based solely on 
RR intervals |33i] . 



D. Systematic error 

There are two intrinsic sources of error in the spectral 
entropy measure related to the phenomenon of spectral 
leakage: that due to the "picket-fence effect" [s^ (where 
frequencies in the power spectrum fall between discrete 
bins) and that due to finite window effects [1£] (where, 
for a given frequency, an integer number of periods does 
not fall into the spectral entropy window). We attempt 
to quantify this error by applying the measure (with pa- 
rameters as per the Data Analysis section) to synthetic 
event series: a periodic series with constant interbeat in- 
terval. For a heart rate range of 50-200 bpm in 1-bpm 
increments we obtain 150 synthetic time series. We find 
the average error in the spectral entropy over the 150 time 
series to be 0.02. The average standard deviation value 
(with variance windows having M equal to 20 spectral 
entropy values) over the 150 time series is 0.011±0.009; 
the average error on these standard deviation values due 
to windowing is 0.0002. 

The presence of some form of error associated with 
finite windows is unavoidable. We have attempted to 
minimize such errors by choosing parameters that achieve 
a balance between usability and error magnitude. There 
is still scope for fine-tuning parameters - in particular, 
trying a variety of window shapes to further reduce the 
affect of spectral leakage. However, we find the general 
results to be robust to a range of window parameters, 
implying any practical effect of windowing errors to be 
minimal when compared to the other issues discussed in 
this section. 



V. FURTHER WORK 

Additional directions for this work include refining and 
extending our cardiac study with a view to clinical im- 
plementation. Furthermore, we suggest that rhythmic 
signals arising from other biological systems may have 
application for the techniques described in this paper. 
An investigation of the optimal windowing parameter set 
would be instructive since our findings suggest the exis- 
tence of physiological thresholds in the spectral entropy 
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level and variance that are applicable to a variety of pa- 
tients. As noted at the end of Sec. IIV C[ one challenge 
would be to investigate and improve the utility of the 
measure (alone or combining methods) when applied to 
patients that demonstrate a mix of different pathologies 
and arrhythmias. Adjusting the spectral entropy win- 
dow to covary with instantaneous heart rate so that a 
always contains ten beats exactly would further reduce 
issues related to variations in the heart rate. Extending 
the algorithm to include other dimensions in the disorder 
map (e.g., heart rate) will likely improve the accuracy of 
results and may increase the number of arrhythmias the 
spectral entropy can distinguish between. 

An accurate automatic detector of atrial fibrillation 
would be clinically useful in monitoring for relapse of 
fibrillation in patients and in assessing the efficacy of an- 
tiarrhythmic drugs [35| . An implementation integrated 
with an ambulatory ECG or heart rate monitor would be 
useful in improving the understanding of arrhythmias on 
time scales longer than that available using conventional 
ECG analysis techniques alone. 

Measures of disorder in the frequency domain have 
practical significance in a range of biological signals. 
For example, the regularity of the background electroen- 
cephalography (EEG is the measurement of electrical ac- 
tivity produced by the brain as recorded from electrodes 
placed on the scalp) alters with developmental and psy- 
chophysiological factors: some mental or motor tasks 
cause localized desynchronization; in addition, the back- 
ground becomes more irregular in some neurological and 
psychiatric disorders ([l0| and references therein). The 
spectral entropy method and the concept of the disorder 
map described in this paper are not cardiac specific: it 
would be instructive to adapt these ideas to other rhyth- 
mic signals where a rapid detection of arrhythmia would 
be informative. 



VI. CONCLUSION 

In this paper we have presented an automatic arrhyth- 
mia detection algorithm that is able to rapidly detect the 



presence of atrial fibrillation using only the time series of 
patients' beats. The algorithm employs a general tech- 
nique for quickly quantifying disorder in high-frequency 
event series: the spectral entropy is a measure of disorder 
applied to the power spectrum of periods of time series 
data. The physiologically motivated use of the spectral 
entropy is shown to distinguish atrial fibrillation and flut- 
ter from other rhythms. For a given set of parameters, we 
are able to determine from a disorder map two threshold 
conditions (based on the level and variance of spectral en- 
tropy values) that enable the detection of fibrillation in a 
variety of patients. We apply the algorithm to the MIT- 
BIH atrial fibrillation database of 25 patients. When the 
algorithm is set to identify abnormal rhythms within 6 
s it agrees with 85.7% of the annotations of professional 
rhythm assessors; for a response time of 30 s this becomes 
89.5%, and with 60 s it is 90.3%. The algorithm provides 
a rapid way to detect fibrillation, demonstrating usable 
response times as low as 6 s and may complement other 
detection techniques. There also exists the potential for 
our spectral entropy and disorder map implementations 
to be adapted for the rapid identification of disorder in 
other rhythmic signals. 
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